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Abstract — The problem of inverting a system in presence of 
a series-defined output is analyzed. Inverse models are derived 
that consist of a set of algebraic equations. The inversion is 
performed explicitly for an output trajectory functional, which 
is a linear combination of some basis functions with arbitrarily 
free coefficients. The observer canonical form is exploited, and 
the input-output representation is solved using a series method. 
It is shown that the only required system characteristic is 
observability, which implies that there is no need for output 
redefinition. An exact inverse model is found for linear systems. 
For general nonlinear systems, a good approximation of the 
inverse model valid on a finite time interval is found. 

I. INTRODUCTION 

Generally, a dynamical system is described by a model 
consisting of a set of differential equations. With knowledge 
of the control input, the states and the output are computed 
by solving these equations. For many controller designs, 
however, a so-called inverse model is required. Here, with 
knowledge of the output, the states and the input are com- 
puted, without solving differential equations. 

An intuitive application of an inverse model is a feed- 
forward controller, where the output reference trajectory is 
given and the corresponding control input is generated. Other 
applications include trajectory tracking control, two-degree- 
of-freedom control, as well as trajectory optimization and 
inversion-based predictive control. 

The methods presented in this paper are applicable if 
the output trajectory is defined as a series functional. In 
many applications, for example in robotics, it is common 
practice to design a continuous trajectory as power series 
through polynomial interpolation between desired setpoints 
[1]. As this is a simple, straightforward and efficient method, 
and as often only the setpoints are of interest but not the 
path itself, this procedure has become popular also for 
feedforward controller trajectories. It is also common in 
trajectory optimization, where the path between the setpoints 
is also considered. Applying the Ritz-Galerkin method (or 
basis function approach), the series definition of the output 
trajectory transforms the optimal control problem to a finite- 
parameter optimization problem [2] [3]. 

Thus, inverse models and series-defined output trajectories 
are very common in the current state-of-the-art as well as in 
research in automatic control. 

This work was supported through the National Research Funds of 
Luxembourg under grant PhD-08-070. 

J-F. Stumper and R. Kennel are with the Institute of Electrical Drive 
Systems and Power Electronics, Department of Electrical Engineering 
and Information Technology, Technische Universitat Miinchen, D-80333 
Munich, Germany, jean-f rancois . stumper@tum.de 



It is noted that the term system inversion is slightly 
different as in the famous works of Devasia [4], where 
numerical integration is involved. The aim here is to provide 
an inverse model in algebraic form. 

Typically, the system is first inverted, then the output is 
defined as series functional [5]. Inversion is then related to 
differential flatness (or feedback linearization), where the 
output is redefined such that the system has full relative 
degree [5]. However, this redefinition, if possible, may be 
undesirable if the output is dependent on uncertain parame- 
ters [6]. Furthermore, tracking of the original output is made 
difficult because of the dynamics between the original and 
the redefined output, which must be respected for reference 
recalculation [7] [8]. In contrast, this method mostly relies 
on observability, and simplifies the implementation as an 
arbitrary output is used. 

For conciseness, only single-input single-output (SISO) 
systems are considered. However, the methods are directly 
extendable to multi-input multi-output systems. 

The paper is organized as follows. In section II, the 
considered problem of inversion with series-defined outputs 
is formulated. In section III, the basic idea of exploiting 
observability is introduced. Section IV shows the method for 
linear systems, and section V treats nonlinear systems. Two 
examples are shown for linear systems and one for nonlinear 
systems. 

II. SYSTEM INVERSION WITH SERIES-DEFINED 
OUTPUTS 

Usually, exact system inversion is achieved through ex- 
ploiting differential flatness, which can be seen as an ex- 
tension to controllability [5]. The result is a differential 
parameterization of the system states and the control input 



x(t)=A(y(t),y(t),...,y^(t)) 
u(t)=B(y(t),y(t),...,/"\t) 



(1) 

(2) 

Knowledge of the output trajectory y(t) leads to the state 
and control input trajectories, without solving differential 
equations. In many applying controllers, subsequently, the 
output trajectory is defined as a series functional 



i=0 



(3) 



to reduce the trajectory problem to a simpler finite-parameter 
problem. The arbitrary free parameters a, could be deter- 
mined by polynomial interpolation between setpoints, or by 
an optimization scheme for trajectory generation. 



In this paper, the two steps of computing the inverse 
and defining the output trajectory are interchanged. Given 
a trajectory defined according to ||3}, the corresponding state 
trajectories x(t) and control input u(t) are searched. Thus, 
the state and control input trajectories are determined as 
algebraic functions of time and the arbitrary parameters a 

x(t)=A(a,t), (4) 
u(t)=B{a,t). (5) 

This is the general form in which inverse models with series- 
defined outputs are applied. The difference in the procedure 
is of interest as the system class is not limited to differentially 
flat or controllable systems, and as the inversion is performed 
for an arbitrary output. 

III. PARAMETERIZATION BASED ON 
OBSERVABILITY 

Definition: Observability 

A system is observable if for any T > it is possible to 
determine the state of the system x(T) through measurements 
of y(t) and u(t) on the interval [0,T]. 

Thus, exploiting the property of observability for system 
inversion seems a good choice. The next two subsections 
present system-theoretic results that, in the following, lead 
to a solution of the inversion problem. 

A. Parameterization of the system states x(t) 

The observer canonical form for a nonlinear SISO system 
[9] is given as 

x\ = x 2 +g\(xi)u 
X2 = x 3 +g2(xi,x 2 )u 

< : (6) 

x n —\ — x n + g n —\ (x\ , . . . )m 
^ x n =F(x) +g n (x u ...,x n )u 

where the output is y = x\. It is noted that, however, the 
observer canonical form can not be found for every nonlinear 
observable system. A parameterization of the system states 
through the output y and input u is found straightforward 
by successive elimination of the states x,-. The resulting 
parameterization is 

X\ =y, 

X2=y-gi(y)u, 

*3 = 1*2- g2(y,X2)u, (7) 
Xn = dx x n— 1 — 8n— 1 ( x l > ■ ■ ■ >X n — l)u. 

Computation of these equations is straightforward. Note that 
the last equation of © is unused. 

A point to consider here are initial conditions. If a system 
is of full relative degree, the initial conditions are regarded on 
the output trajectory only. If not, as in the examples in this 
paper, some conditions involve the control input trajectory 
u(t) too. 



B. Parameterization of the control input u(t) 

A further result on observability is that for every observ- 
able system, an input-output representation can be found 
[10]. Input-output representations are higher-order differen- 
tial equations in the inputs and outputs, and are equivalent 
to state-space descriptions. The input-output equation takes 
the form 

q(y,..., yW «,..., «(»)) =0, (8) 

in which the output appears with its n-th and the input 
with its m-th derivative. To determine the input u(t) for a 
given output trajectory, ^ is solved. It is noted that this is 
equivalent to solving the last line of (0. 

Thus, for a full parameterization of the system variables, 
the solution of one differential equation is required. If the 
output is given as (O, the solution of the differential equation 
is possible in a straightforward manner with the widely 
known series methods [11]. The solutions are presented in 
the next two sections. 

IV. RESULTS FOR LINEAR SYSTEMS 

This section demonstrates that an exact inverse model is 
easy to compute for an arbitrary output of a linear observable 
system. 

A. Exact inversion with a power-series output 

The state trajectories x(t) are determined by successive 
elimination of (0. As the observer canonical form is linear, 
the functions gj(-) are constants and the function F(x) = q T x 
is linear. The result is 

x\ =y, 
X2=y-giu, 

X3=y-giu-g 2 u, ( 9 ) 

x n =y("- l '>-g l u(>-V-...-g n _ l u. 

Furthermore, the input-output representation is a linear equa- 
tion, as follows from the last equation of (0: 

x n = q T x + g„u, (10) 

and consequently, 

q T x = yW-g l u("-V-...-g n _ l u-g„u 1 (11) 

which becomes a linear higher-order differential equation in 
u and y with ®. For conciseness, it is rewritten as 

n m 

EAw.v' (12) 
The output trajectory y(t) is defined as a power series 

N 

y{t) = Y J a i t i , te [t ,t f ], (13) 
;=o 

meaning in © is set *¥i(t) — t'. Power series simplify the 
procedure, as ^ ^ (f ) = i (f ) and 'Po = 1 ■ In the examples, 
it will be shown that this choice is useful for observable and 



controllable systems, but that for non-controllable systems, 
other basis functions shall be used. The derivatives are 

y {,) (t) = j^-aj+itJ, t g [t ,t f ). (14) 

The applied method is known as method of undetermined 
coefficients [11]. Explicit computation of the initial condi- 
tions u(to),..,u( m ~ l '(to) is avoided. The control input is a 
power series with the undetermined coefficients /}, 



«(0 = E# f '' f e M/]- 

i=0 

Both trajectories are inserted in (fT2l to obtain 



(15) 



A?—! 



!=1 /=QU v- 



7+' 



(16) 



The undetermined coefficients j3,- are then found comparing 
the terms in equal power of t. To find the inverse model, 
the parameters /3, are eliminated and replaced in (Tl3T > by the 
linear functions of 05. 

Exactly n of the coefficients OS* are fixed by the initial 
conditions x(t$). Evaluation of (0 leads to the respective 
equations. It is however noted that the information included 
in (0 is also included in the input-output representation, thus 
the initial value problem could be transformed such that some 
of the coefficients j3; are fixed. 

The main result is that the differential equations are ex- 
actly solved. The inversion is equivalent to an exact inversion 
followed by an output trajectory approximation ([3). 

B. Example of a linear controllable system: Buck converter 

The average model of the buck converter, or synchronous 
step-down converter, is given by the linear lifferential equa- 
tions in the capacitor voltage u c and the inductor current // 



U dt Mc — {R c +R) Uc 



*4« = 



(R C +R) c 



R ■ 

(Rc+R) 1 ' 

Rl 



RRc \ 
(Rc+R) ) 



h + Uitfii 



, (17) 



where the control input u is the duty cycle. The output is the 
voltage y = Uqut 



Uqut = 



R 



( ggg 

R c + R) ut \Rc + R 



(18) 



It is noted that u c is the controller canonical form output, not 
Uout- For conciseness, the parameters are set to Um = 12V, 
C = 200AiF, L = lOOjiiH, R = 2Q., R c = 85mi2 and R L = 
280mi2. The states are eliminated via the observer canonical 
form, yielding 



u c = 0.98y - 18.6y+ 181950m, 
U =0.707y + 219y- 2143200m. 

The input-output representation therefore is 

y + 0.006 y = 575m + 9780m. 



(19) 
(20) 



(21) 



The output trajectory and the control input are defined as 
degree 3 power series 

3 3 

y(f) = »W = IM fg[0,T]. (22) 

1=0 i=0 

To determine the coefficients j3,-, the linear system of 
equations when inserting (1221 in ( f2TT > is solved. Comparing 
the terms in equal power of f, the undetermined coefficients 
are found as 



/3 = 0.00001043ai + 0.003123a 2 -0.1593a 3 , 
j3i =0.00002086a 2 + 0.009369a 3 , 
j3 2 = 0.00003130a 3 , 
J8 3 =0. 



(23) 



Setting to ~0, the initial conditions are m c (0) and i c (Q). 
Evaluating (fT9l and (l20l fixes the first two parameters of the 
output trajectory ab and a\ as 

ao = 0.9615« c (0) + 0.08167z/(0) + 0.3065O2- 15.64a 3j 

(24) 

Cfi = -0.003457M c (0) + 0.004792i / (0) + 34.0402- 1737a 3 . 

(25) 

The coefficients a 2 an d 0^3 are tnus f ree an d can be use d 
in a trajectory generation scheme, for example reference tra- 
jectory interpolation. For an output trajectory approximated 
according d22L the inverse model is found in terms of the 
arbitrary parameters of the output trajectory a,-. The states 
are 

u c = (0.98oo - 16.70a! +568.2a 2 - 28984a 3 ) 
+ (0.98ai - 33.40a 2 + 1704a 3 )f 



(0.98a 2 -50.10a 3 )f 2 + 0.98a 3 f 3 



(26) 



ii = (0.707 ao + 196ai -6693a 2 + 34141 la 3 ) 
+ (0.707 a, + 393.2a 2 - 20079a 3 )f 
+ (0.707 a 2 + 589.9a 3 )f 2 + 0.707 a 3 f 3 , (27) 

and the control input is 

m = (1043cei + 31 .23a 2 - 0.1593a 3 ) + (2086a 2 + 93.69a 3 )f 
+ 3130a 3 f 2 . (28) 

C. Example of a linear uncontrollable system 

To demonstrate that the proposed inversion scheme is only 
based on the observability of a system, the method is now 
applied to an uncontrollable and thus also non-flat system. 
The constructed example is defined as 



X\ — X 2 + M 

±2 = -2x 2 



(29) 



with the output y =x\. The system is in observer canonical 
form, and has an uncontrollable subsystem x 2 . The input- 
output representation is given as 



y + 2y = m + 2m, 



(30) 



where a priori, no pole-zero cancellation happens. If power 
series are applied, some terms will cancel in the input- 
output representation, same as in a transfer function in 



Laplace domain. The initial conditions would be restricted to 
^2(0) = 0. To avoid this, the basis functions must be chosen 
such that iWi(f) ^ For instance, possible choices 

are exponential functions ^P, = e~" or trigonometric series 
Wi = sin(/f). Also, a combination of such functions with 
power series is interesting, as the uncontrollable subsystem 
and the input-output behavior can be designed independently. 

The output trajectory and the control input are defined as 
the degree 3 exponential series 

3 3 

= £«,•*-''', 11(f) = te[0,T}. (31) 

1=0 i=0 

Placing this definition into the input-output description, and 
comparing the terms in the respective exponent, three unde- 
termined coefficients are found, 

j8o = 0, (32) 
J8i = -ai, (33) 
ft = -3a 3 . (34) 

The initial conditions define the first parameter of the output 
trajectory, as well as the missing parameter of the input 
trajectory 

Ob = jci (0) — oei - a%- 0C3, (35) 
p 2 = -x 2 (0)-2a 2 . (36) 

The inverse model for the given output trajectory is thus 
found as 

xi(t) = ao + aie~' + a 2 e~ 2 ' + a 3 e~ 3r , (37) 
x 2 (t) = -2a 2 e- 2 ' - p 2 e- 2 > = x 2 (0)e- 2 ', (38) 
u(t) = -a x e~ x - (x 2 (0)+2a 2 )e- 2t - 3a 3 e~ 3 ' ■ (39) 

A posteriori, the restriction j3,- = could be removed, for 
instance by adding a term a^t to y(t). This result shows 
that the choice of basis functions VPj(f) is important. Which 
type of functions is advantageous is depending on the system 
characteristics as well as on the application. 

V. RESULTS FOR NONLINEAR SYSTEMS 

Successive elimination of the observer canonical form (O 
yields the parameterization of the state trajectories x(t) 0. 

For a full parameterization, the control input u(t) is de- 
termined by solving the input-output representation ([8]). The 
nonlinearities are inherited therein. Even though an analytic 
solution of this differential equation exists in many cases, 
this is typically an implicit solution in terms of elementary 
functions, which may be less useful than an explicit series 
solution [12]. 

To solve the differental equation, the method of undeter- 
mined coefficients is extended to nonlinear systems [11]. An 
alternative would be the method of successive differentiation, 
which is however only applicable for systems with a well- 
defined relative degree. This method also works with singu- 
larities as «W is not eliminated. Singularities are common 
in non feedback-linearizable (or non-flat) systems. 



With (0), equation ([S]) becomes 

f(a,t,u,..,u^) = 0. (40) 

Assuming that the input u is defined as a series of same type 
as the output y(t) in (fJJ 

N 

«(0 = te[0,r], (41) 

!=0 

the input-output representation is further simplified to 

/(a,j3,f) = 0. (42) 

The function / is supposed to be analytic, meaning, / 
has a Taylor series expansion in powers of f, valid in a 
neighborhood of to- 

As there are N +1 undetermined coefficients j3,-, N +1 
independent equations must be evaluated. In a nonlinear 
system, for instance comparing terms in equal power in 
t when power series are used, as in the previous section, 
is generally not possible as the degree of d42l may be 
higher than N+l. However, a similar result is obtained by 
evaluating the equations 

f®{a,p,t) = if(a,li,t) = 0, i = 0..N-l. (43) 
dt 

Some coefficients j3; may be well defined by the initial 
conditions. For the remaining R < N+ 1 coefficients, a Taylor 
series solution of d42l around a point t = to is searched. The 
remaining R undetermined coefficients j3; are obtained by 
solving the system of algebraic equations 

/(Cf,j3,fo) = 0, 

. (44) 

f( R -V(a,p,t Q ) = 0. 

which is equivalent to solving the input-output representation 
by ignoring the higher-order terms. This method is compa- 
rable to Taylor series approximation of the input trajectory. 
As this is a system of nonlinear equations, finding a solution 
may be difficult. Thus it cannot be guaranteed whether the 
method is applicable for a given system. 

The inverse model is exact at t = to, but generally has 
an error for t ^ to- The error is bounded in a finite interval 
around to and is dependent on the order of the series N. Also, 
the choice of basis functions *Pf(f) is relevant. Thus, only an 
inverse model of the type (01 (0 for 'short' time intervals 
around t = to can generally be found. 

An open question is the convergence of the input trajectory 
to the exact input, which is of interest to determine the 
required number of parameters /3,. The number of parameters 
for a finite interval may grow exponentially with the interval 
length. 



A. Example for approximate inversion: Van de Vusse reactor 

The Van de Vusse type reactor is presented in [13] as a 
benchmark example, mainly for operating point changes of 
nonlinear systems. The state-space model is 



x\ = —k\X\ — k^x\ + u(Caq— Xi) 



(45) 



X2 = k\X\ — k2X2 — UX2 

where X2 represents the product concentration in the output 
stream (mol/L), and x\ the reactant concentration inside the 
reactor (mol/L). The control input u is the dilution rate in 
the input flow relative to the reactor volume (1/h), and Cao is 
the reactant concentration in the input system (mol/L). The 
output is the product concentration 



y=*2, 



(46) 



and the system is observable. Even though the system is 
differentially fiat [14], a flatness-based design requires output 
redefinition to a full relative degree output. As X2 has relative 
degree one, there are internal dynamics associated to this 
output. 

For conciseness, the parameters of the system are fixed to 
kt =50i, k 2 = 100 J, h = 10^, C A0 = 10^ and V = IL. 
The states are eliminated as 



*2=y, 

x\ = — y + 2vH yu, 

50 50 

and the input-output-representation is 



(47) 
(48) 



50y + 75003) + 250000y 
+u(100y 



I0f + 2000xy +100000/ 
7500y + 20yy + 2000y 2 - 250000) 
+u 2 (l0y 2 + 50y) + ii(50y) =0. 



(49) 



The input-output representation is a nonlinear function as u 
is involved as square. The output trajectory and the control 
input are defined as power series 



y 



N 
;=() 



(50) 



= £M te[t ,tf], 

i=0 

with an output trajectory degree N = 3 and with input 
trajectory degree N' to be determined later. The system of 
equations ( f4~4~T > can now be established. 

Setting to = 0, the initial conditions of the system are x\ (0) 
and #2(0). Comparing to d47l i and j48l . the first parameter 
of the output trajectory and of the input trajectory are fixed: 



«0 : 

ft 



: *2(0), 

:50 ^(o) 



a_ioo. 



(51) 
(52) 



Oo Oo 

1 ) Solvability issue: To analyze the issue of solvability, it 
is seen that the undetermined coefficients J3, are involved as 
shown in table U for to = 0. 

As the remaining coefficients Pi.-Pn are sequentially in- 
volved as first-order (j3,- in /'')(•)), they can be determined by 
successive elimination of the equations d44l >. The resulting 
parameterization is available in the desired format u(t) = 
B{a,t) , leading to x(t) =A(a,t) with (08]l. The results are 
not printed due to space constraints. 



TABLE I 

Degree of the parameters in the respective equations 





ft 


ft 


ft 


ft 


ft 


f{a,P,to) 


2 


l 











/ (1) («,/Uo) 


2 


1 


1 








/ (2) («,/Mo) 


2 


2 


1 


1 





/ (3) («,/3.'o) 


2 


1 


1 


1 


1 



2) Quality of the inverse model: The parameterization 
does not solve the input-output representation for all times 
t, but merely for t = to. This is related to the nonlinearity 
of the equation, as in (49), the order of the input-output 
equation is seen to be N' 2 x A^ 2 because of the term u 2 y 2 . To 
solve the system, the higher-order terms are ignored. Thus, 
the found inverse model is not exact and the error or offset 
is proportional to the order of the Taylor polynomial N' in 
d50l . If more information is desired, N' must be increased. 

A steady-state setpoint cange from X2 = 0.9 to X2 = 1 • 1 in 
T = 1 h is considered. This fully determines the parameters 
OC,-, and the initial conditions define j3o- Furthermore, the 
system is in minimum-phase operation. The output recon- 
structed from the control input u(t) by Euler integration 
is compared to the original output trajectory. Results for 
different orders N' are displayed in Fig. [T] namely N' = 3, 
N' = 5 and N' = 9. In the results, x\ (?) is always good 
as the impact of u(t) is smaller than the impact of the 
output trajectory. The input, however, is converging slowly 
to the exact inversion input for an increasing order N'. For 
a quantitative comparison, the integrated trajectory offset 



E = 



y(t)-y(t) I dt 



(53) 



is shown in table [TT] Although it seems that the solution 
is exact for N' = 9, it is actually not, as there are higher- 
order terms in (50) as mentioned before. The found inverse 
models are usable only on a give time interval t 6 [fo,f/]. Still, 
the results for N' — 9 are exact enough for an application, 
as other error sources, such as parameter uncertainties, may 
have more impact as the error of the inverse model. 

TABLE II 

Trajectory error for the inverse models 



degree N' 


3 


5 


9 


error E 


0.0964 


0.0050 


0.00006 



VI. CONCLUSIONS AND FUTURE WORKS 

A. Conclusions 

In this paper, the problem of system inversion in the 
presence of a series-defined output trajectory is analyzed. 
Overall, the problem of system inversion is greatly simplified 
if the output is considered as a functional, which is a 
linear combination of basis functions with free coefficients. 
Known results on observability are used to establish ordinary 
differential equations. Subsequently, these are solved using 
series solutions and the method of undetermined coefficients 



1 1 1 1 1 1 1 original trajectory 




time t [h] 

Fig. 1. Van de Vusse type reactor: Results on inversion. Original product 
trajectory X2, recomputed product trajectories from u and {45) (Euler 
integration), reactant x\ and input dilution u computed by inversion, with 
N' = 3,5,9. 

[11]. An (algebraic) inverse model is found for an arbitrary 
given output, it consists of a set of equations for the system 
states and control input depending on time as well as on the 
(free) output trajectory coefficients. 

An exact inverse model is found for linear systems. For 
general nonlinear systems, assuming that a system of given 
nonlinear equations can be solved, a good approximate 
inverse model valid on a finite time interval is found. 

The presented methods, compared to the existing 
differential-algebraic inversion schemes, have the advantage 
that the system is not required to be differentially flat (lin- 
earizable by feedback), and that the output is not redefined 
to the controller canonical form output [5]. Compared to 
the existing numerical schemes [4], no differential equations 
are numerically integrated, which is an advantage if the 
system is of a high order or the trajectory horizon is high. A 
comparison to numerical methods is however not very useful, 
as the applying controller designs are different. 

B. Future Works 

For many applications, especially of the class of non-flat 
systems, it is hard to obtain an inverse model or to design 
a feedforward controller [5]. If such a system is observable, 
it should be analyzed whether the method yields a suitable 
solution valid on a practicable time interval. Furthermore, 
in terms of robustness, it should be analyzed whether an 
approximate model is better than one with a redefined output. 
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